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^) • The infinitesimal unitary transformation, introduced recently by F. Wegner, to bring 

the Hamiltonian to diagonal (or band diagonal) form, is applied to the Hamiltonian theory 
as an exact renormalization scheme. We consider QED on the light front to illustrate the 
method. The low-energy generated interaction, induced in the renormalized Hamiltonian 
to the order a, is shown to be negative to insure, together with instantaneous term and 
perturbative photon exchange, the bound states for positronium. It is possible to perform 
the complete elimination of the ee7~vertex in the instant form frame; this gives rise to the 
cutoff independent ee-interaction governed by generated and instantaneous terms. The 
well known result for the singlet-triplet splitting |a 2 Ryd is recovered in the nonrelativistic 
limit as long as A <K m. 

We examine the mass and wave function renormalization. The ultraviolet divergencies, 
associated with a large transverse momentum, are regularized by the regulator arising from 
the unitary transformation. The severe infrared divergencies are removed if all diagrams 
to the second order, arising from flow equations method and normal-ordering Hamiltonian, 
are taken into account. The electron (photon) mass in the renormalized Hamiltonian vary 
with UV cutoff in accordance with 1-loop renormalization group equations. This indicates 
an intimate connection between Wilson's renormalization and the flow equation method. 

The advantages of the method in comparison with the naive renormalisation group 
approach are discussed. 
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1 Introduction 

We apply the method formulated by F. Wegner |l| of unitary transformation, to bring the 
Hamiltonian to diagonal (or near diagonal) form, to the problem of renormalization. We start 
with the regularized bare cutoff Hamiltonian, with the energy widths (the energy differences 
between the free states in a matrix element) restricted to be below the bare UV cutoff A. As 
renormalization step we perform the unitary transformation that removes the energy widths 
in the final Hamiltonian, to be below the final cutoff A < A, resulting in the band-diagonal 
structure for the renormalized Hamiltonian. This corresponds to the Wilson renormalization 
procedure in the Lagrange theory of integrating out the fields of higher energies. 
This unitary transformation can be written 

U(X, A) H B (A, e A , m A , g(e A , m A )) U + (\, A) = H B (X, e x , m x , g(e x , m x )) , (1) 

where e x ,m x denote the set of finite number of independent coupling constants, in the case of 
QED - the ee7-coupling constant and the fermion mass. Instead we generate in the renormalized 
Hamiltonian new effective interactions, denoted as g(e x m x ), and corresponding to the operators 
of higher dimensions. They must be added to the bare cutoff theory as irrelevant operators, 
namely 

g(e A ,m A ) = , A -*■ oo , (2) 

being manifest at the low-energy scale. The number of independent coupling constants of the 
theory is preserved, i.e. g — > when e — > 0, m — > 0. 

By definition the renormalized Hamiltonian, see , in the literature also called the effective 
Hamiltonian, is 

H R = lim H B . (3) 

A^oo 

For renormalizable theories there arise only a finite number of infinities, associated with 
the limit A — > oo. They can be eliminated order by order in the running coupling e x by an 
appropriate choice of cutoff-dependent counterterms. The counterterms are chosen from the 
constraint on the renormalized Hamiltonian to run coherently with the cutoff. 

This renormalization scheme, called similarity Hamiltonian approach, was presented by 
S. Glazek and K. Wilson || as an alternative to the traditional Lagrange method. 

The similarity unitary transformation was used by these authers in light-front field theory 
to obtain the band-diagonal Hamiltonian with finite cutoff A. For these purposes we exploit 
in this work the flow equations, formulated by Wegner []]], that perform infinitesimal unitary 
transformations in a continuous way. 

We consider in this work two aspects of renormalized Hamiltonian. The first point is related 
with renormalization group analysis. The second one is connected with the low energy sector 
of the theory and bound state calculations. 

In the next section we give the key ingredients of the flow equation method and discuss and 
its application to renormalization group analysis in Hamiltonian theory. 

2 Flow equations 

The flow equation is written 

dJ ^=W),H{l)] (4) 



or for the matrix elements 

^ = r hd (E j -E i ) + [r,,H I ) ij , (5) 

with the Hamiltonian H = H + Hj devided into a free and an interacting part. Here rj(l) is 
the generator of the unitary transformation, that determines in the leading order the evolution 
of the interacting part Hj^, i ^ k or, vice verca, can be regarded as function of Hj, i.e. rj(Hi). 
The flow parameter I changes from l^ = as A — > oo, that corresponds to the initial canonical 
Hamiltonian H can , to some finite value l\, where the Hamiltonian has band diagonal form. 
The value l\ = oo corresponds to the diagonal (or block- diagonal, in the case when exact 
diagonalization is impossible) form of the Hamiltonian. The dimension of the flow parameter I 
is 1/ (energy) 2 . 

In the problem of renormalization the physical sense of the UV cutoff can be assigned to 
the flow parameter I, explicitly it satisfies 

h = y 2 , (6) 

where A is the UV cutoff. The unitary transformation removes the matrix elements between the 
free states with energy differences between the bare cutoff, A, and the final cutoff, A (A < A), 
resulting in a band-diagonal form for H\. 

We continue with the method of flow equations itself. The generator of the transformation, 
according to 0, is chosen 

Vij = [Ho, H}ij = (Ei - Ej)Hi tij (7) 

giving 

*!k = _(£. _ E 3 f E lM + fo, Hfa . (8) 

Generaly the interacting part is associated with the coupling constant g. We proceed further 
in the frame of perturbative theory (PT) with respect to the interaction Hi, or therefore 0(g). 
From the equations above follows, that the generator rj to the leading order 0(Hj) generates the 
effective interaction to the order O(Hf). The new term must be added to the initial Hamiltonian 
as irrelevant operator at /a — ► (i.e. vanishing at I = 0, H gen (l = 0) = 0) but being nonzero 
for any finite value of l\. To make the effective interaction band-diagonal we need to introduce 
the generator to the next order 0(H]), which in its turn generates the effective interaction to 
the third order 0{Hj). We truncate this series assuming the coupling constant to be small. 
Therefore for finite value of l\ one has 

H{1) = H (l)+H<U{l)+HW(l) + ... 
V {1) = r,U(l)+ V W(l) + ... , (9) 

where the upper index denotes the order of PT, and H^(l) = Hj(l), with liniA^oo Hj(l\ = 0), 
being the interaction present in the initial canonical Hamiltonian. 
To the leading order of perturbative theory (PT) 

d -^ = -{E l -E j y H ^. (10) 

Neglecting the dependence of the energies Ei on the flow parameter (being of higher order), we 
obtain the solution in the form 



H$%) = H It «{l x ) = iWU)^7TT , (HJ 



where we have intoduced 



f l3 (h) = exp {-l x (Ei - E 3 ) 2 } = exp j 



Ei~-E]j x 
A"" 



(12) 



The initial condition is defined at the bare cutoff A — > oo, namely lim^^ Hj(lh) corresponds 
in the case of QED to the bare coupling constant e . 

The structure of the expression for Hj(l\) is transparent. By the unitary transformation 
performed with the leading order generator 77, we succeded to eliminate far-off-diagonal elements 
from the Hamiltonian matrix, so that only matrix elements with energy differences \E{ — Ej\ < 
4= = A are present. 

To the next to leading order 0(H 2 ) 



dH 



(2) 



V_ - _(W. _ EM2rr(2) 



-(E.-E^H^ + i^^H^. (13) 



dl 

We introduce H\»- (l\) = fij{l\)H\ - (l\), then the flow equation is rewritten (again neglecting 
the cutoff dependence of the energies Ei) 



dHi; 



(2) , 

This gives rise to the following solution 

'(/a) = H\f(l K ) + f X dl'h^\H^] l3 . (15) 

Jl\ Jij 



hP> 



We are not going to give further details. But it is worth to mention here, that in the case 
of QED the second order term H ^ contributes (in different sectors) to the new generated 
interaction H^ gen ^ and to the electron (photon) self energy terms. This means that beginning 
from the second order the values H^j(lx), defined as H;" (l\) = fij(l\) H^ (l\), contain two 
terms, the interaction part and counterterms. Generally it can be written 

Hto{l x ) = Hl n \l x )+XW{l x ). (16) 

The counterterms are determined from the coupling coherence condition, 0. For the explicit 
treatment of the second order for QED, see later. 

Generally, the system of self-consistent flow equations for the Hamiltonian H(l) and the 
generator of the unitary transformation 77 (Z) can be written in matrix form as follows 

dH i:j (l) _ ^ trl t dlnfij 

(17) 

where the function fy determines how fast the non-diagonal part of the Hamiltonian matrix 
with \E{ — Ej\>\ vanishes with I, /^ = exp(—l(Ei — Ej) 2 ) Note, that this form of the function 
fij corresponds to the choice 77 = [Hq, H] for the generator. Other choices for the function fij(l) 
are possible (see Appendix A). 

The structure of the flow equations, eq. (|T7|), is transparent. The first term in eq. ([17]) 
is responsible for the renormalization of physical values (coupling constants, masses, wave 




functions) and also for the structure of the new generated interactions to the higher order (with 
respect to the canonical interaction) being then present in the initial Hamiltonian; the second 
term in eq. ( |T/j) insures band diagonal structure of the renormalized Hamiltonian. As will be 
shown later, in the case of QED the physical masses run to the order g 2 , while the coupling 
constant starts to run to the order g 3 with the right renormalization group coefficients, i.e. with 
those obtained in standard perturbative theory. 

The second equation for rj is chosen in the form that the terms with small energy denom- 
inators (\Ei — Ej\ — > 0) effectively do not contribute in the renormalized Hamiltonian. This 
advantage, as compared with naive perturbative theory, enables to consider Hamiltonians with 
continuum spectrum. We mention that the problem of small energy denominators is solved 
also in other schemes (Appendix A). 

3 Renormalized Hamiltonian Hr to the second order 

3.1 Canonical light-front QED 3+ i Hamiltonian 

We start with the canonical light-front QED Hamiltonian H can , devided into free and interacting 
parts 

P- = H can = J dx-d 2 x L {Hv + Hi) . (18) 

In light-front gauge A + = A + A 3 = 0, the constrained degrees of freedom A~ and ip- (ip = 
ip + + tp_,ip± = A±ip) can be removed explicitly; this gives the light-front gauge Hamiltonian 
defined through the independent physical fields A± and ip + only || 



H = -(d i A j )(d i A j )+t 



-d\ + m 2 " 
id+ 



£ 



'I'l iLecy "i '^ecyy i" 'f- 



and 



H 



eej 



e? H 



-2§-A-) + *-A^-^ + °-^*-Ai 



n 



ee-yy 



-ie 



H e 



2e 2 



eo-A^-{a-AH) 



(19) 
(20) 

(21) 
(22) 



(23) 

where {a 1 } are the standard 2x2 Pauli matrices, and d + = 2d- = 2^|^. We have used 
the two-component representation for fermion fields introduced by Zhang and Harindranath 
H ip + = (q). To simplify the calculations we rewrite all interactions through creation and 
annihilation operators. This turns out to be useful in the flow equations formalism, j3]. 

Following standard quantum field theory procedure we use the momentum-space represen- 
tation for the field operators, M and H, 



A\x) 



Ex,/ 



E 

A 



dp + d 2 p J 
2(2tt) 3 
dq + d 2 q ± 0(q 

2(2tt) 3 



e(p + ){b p ^ x + d p , s e i n 



(s\a q , x e'^ + h.c. 



(24) 



where spinors are xXii = (1, 0), X-\ii tr = (0,1), with s = — s and polarization vectors e\ = 
5=(l,i), e l _ x = 75(1, —i)', the integration running over the p + > only these states, that are 
allowed the light-front theory. 

The corresponding ( ant i) commutation relations are 

\"p,s) bpigtf = \d PiS , dp/gij = o P: pid ss i 

[°«,A, oJ,A'] = Sq,q'S\,\< , (25) 

where 

S PjP > = 2(2irf 8(p+ - p'^ 2 \p ± - p' L ) . (26) 

The light-front vacuum has trivial structure for both boson and fermion sectors, namely a q \0 >= 
0; b p \0 >= 0, simpifying the analitical calculations. The normalization of states is according to 

< Pi, si\p 2 , s 2 >= S PUP2 8 SUS2 , (27) 

where fr^jO >= \p, s >. 

Making use of the field representation eq. (^), we have the following Fourier transformed 
for 
the free Hamiltonian 

h« = E / d i^»(p + f-^-(KA.* + <A.) + E / %£f g ('^<y^ • < 28 ) 

the leading order 0(e) ee7-coupling 

Hee-y = Z^ / [9p 1 p 2 qV') £ \Q'q + 9pip 2 qV') £ \ &-q\ ("p 2 "Pi + % -pi + °^-P2"Pi + "-P2"-pJ 
As lS2 J P1P2" 

x x^rKpi,P2,-g)x S i^,p 2 -pi , (29) 

where 

Ti{Pi,P2,q) = 2-r + <r-o + , 30 

q+ pj pj 

where the mass is 1-dependent. Further we have for the instantaneous interactions of the 
order 0(e 2 ) 

HTeil = E / <™ 4 (0(^3 + ^ P 3)(^ 4 + rf- P J(V + ^J(^ + ^ P2 ) 

SlS 2 S 3 S 4 J P1P2P3P4 

4 

X Xs 3 Xs 4 7^j: +^2XsiXs 2 0p3+P4,Pl+P2 (."-U 

{Pi ~ P3 ) 

and 
KZ = E / 9 e pJ 2 f qi q 2 (l)^xV< + e Xl ~a- qi )(e{ 2 ~a q2 + e% l a±J(b+ + L P2 )(b Pl + dt Pl ) 

^Xs 2 ~T^ i + \Xs^p\+q 2 ,q\+p 2 , (."^J 

here 

9(-g+) 



3 = n 8 -^-l 
V9+ 



O— o — «— o.A" 

V-9" 1 



& P = b PiS 9(p + ), dp = d PtS 0(p + ) , (33) 



and the 5 stands for the function defined in eq. fl26|), the short notation for the integral we 

understand as 

dp + d 2 p 1 - 

, , 2(2tt) 3 ' l } 

In the formulas above we write explicitly the momentum dependence of the coupling con- 
stants as long as I ^ 0. The initial conditions for the couplings are defined at the value of the 
bare cutoff A — > oo, namely 

lim g ee \l k ) = e (35) 

A— >00 



and for both instantaneous interaction couplings 

lim g mst (l A ) = e 2 ; (36) 

A^oo 

these correspond to the couplings of the canonical theory. 

3.2 The flow equations in |ee>-sector 

3.2.1 Generated interaction 

Following the procedure outlined in the second section, the leading order generator of the 
unitary transformation is 



x xt 2 r i(PuP2, -q)xs 1 s q , P2 - P1 , (37) 

where Pi and pj stand for the set of initial and final momenta, respectively, and 



VpiPf{'') — ^PiPfdpiPf — a ' ,n ' y 



Further we calculate the bound states of positronium. In what follows consider in \ee > 
sector 
the generated interaction to the first nonvanishing order 

-"eeee ~ Z_^ j PiPf \ ' V3 PA P2^PiXs 3 Xs 4 Xs 2 Xsi0pi+p 2 ,p 3 +p4 ) V*™) 

S1 S 2 S 3 S4 J P1P2P3P4 

with the initial condition liniA^oc V^{1\) = 0, 
and the instantaneous interaction 

-"eeee = 2-^i / PiPf w"p3^P4 < ^P2"PlXs 3 Xs 4 X52Xsi ( 'pi+P2,:P3+P4 ) l^U) 

S1S 2 S 3 S4 J P1P2P3P4 

where 

VZl(l)=9 P X(l)j-T A - f y(l)- (41) 

The ordering of the field operators in both interactions is important because it defines the kind 
of interaction, attractive or repulsive; the ordering given satisfies the standard Feynmann rule 
prescription in the |ee > sector. 

We neglect the I dependence of momenta in the interaction, which enables us to write the 
flow equations for the corresponding couplings. 



The flow equations to the order 0(e 2 ) 
dg P iPf{l) 



dl 
dl 

dv P T f ( l ) 

dl 



where 



PiPfSpiPf^) 

<[ V ^(l),H eey ]> ]e - e> -Al iPf V p %(l) , 



A PrPf =J2Pi ~Y,Pf 



(42) 



and the light-front fermion energie is p 



p ±2 +m 2 



the photon one q 



2Z 



(43) 
The matrix 



element < [q^(l), H eei ) >\ eS> is understood as the corresponding commutator between the free 
electron-positron states, namely < P2,Sz ) PaSa\---\piSi 1 P2S2 >■ 

Neglecting the dependence of the momenta (the fermion mass) on the flow parameter /, the 
solution reads 



9piPf\< j ) 


= f PiPf ■ e + 0(e 3 ) 


C1(0 - 


= W^ + 0(e 4 ) 


V9 en (l) = 


r l l 



(44) 



JpiPf (.' 



•\ee> 



where the subscript indicates, that the commutator is considered in the electron-positron sector. 
As was discussed in the previous section, the function f PiP , eq. (|12D insures band-diagonal 
structure for the renormalized interaction. 

The matrix element of the commutator [i]^\ H eei ] in the exchange and annihilation channels 
is (Appendix B) 



< [r] {1 \H eej ] > /8 P1 . 



-P2,P3+P4 



^2ii f p +_p+\\Vpi,P39p4,P2 I" Vp4,P29pi,P3J > 



M-2U t p + +p +\ \ r lpi,-P29p4,-P3 I" Vp4,~P39pi, -P2 ) 5 



(45) 



where 



Op 



1-P2 



e - 



1 dl 



P1,P2 



PlP2 "" 



ft»l 



P2 



eo •/. 



P1.P2 



(46) 



P1.P2 



Pi — £>2 — (pi ~~ P2) ■ In what follows we drop the index at the coupling. The 



and A 

matrix elements M2U between the corresponding spinors in both channels 



M. 



(ex) 



M. 



2y 

(an) 



[xt 3 r i(P^P3,Pl -P3)Xei] [Xfji(-P4,-P2,-(P1 -P3})X 



*4J 



(47) 



2iJ 



[xt 3 r i(-P^P^ ~(Pl +P2))Xs 4 ] b&FliPl' -P2,Pl+P2)Xsi] 



determine the spin structure of the effective interaction. It is calculated explicitly in light-front 
frame in Appendix B. 



The form of the second order renormalized interactions eq. (f|§), expressed through the /- 
function (that defines the behaviour of the first order coupling), is universal for the different 
renormalization schemes (see Appendix). (This is true up to some factor in the integral for the 
generated interaction). Specifying the / function we obtain the explicit form of the renormalized 
interactions for the different unitary transformations. In Appendix C we compare the results 
for the second order generated interaction in two renormalization schemes. 

Here we choose the /-function as 

f Pi , Pf = eM-lK, Pf )- ( 48 ) 

Then, neglecting again the dependence of the momenta on the flow parameter I, we have 
for the generated interaction in both channels 



llJpitPZV 1 J/P4,P2V' ) 



1 fi 

V ex (I) = -e M 2ii — p -T-(A PuP3 + A P4tP2 )f PiPf - dl 

(Pi ~~ Ps ) J0 JPiPf( L ) 

(Pi + P2 ) J0 



(49) 



JPiPfK'' ) 



This gives rise to 



9en{) 2n (pt- P t)2W,Ps a P4 , p J y 1 ) 



y(an)(j\ _ p2 M (an) 1 1 ( \ l \ /, 

gen {) ~ 2U (Pt+Pt) 2 U P1 ,- P2 V-J [ 



(50) 



\i 



e -/-2A Pli _p 2 Ap 4i _ P3 \ _ g l/ ^PiPf 



where 

^PiP/ —Pi +^2 _ ^3 _ ^4 = ^Pl,P3 ~~ ^P4,P2 = ^Pl,-P2 _ ^P4,-P3 (5*-) 

due to momentum conservation in '+' and 'transversal' directions. 

For energy conserving processes, i.e. when A PtPf = 0, the unity in eq. ([51]) gives the result 
of the standard renormalization procedure. This interaction containes divergencies in the form 
of small energy denominator, as is general for perturbative approach. This problem is cured 
in the method of flow equations (and also by similarity transformations) by the proper choice 
of the generator rj (Appendix C). The divergencies in eq. (|5T|) effectively are cancelled by the 
exponential factor in the bracket (1 — exp) as long as the leading order ee7 coupling is not 
completely eliminated (i.e. for the finite cutoff l\). 

We rewrite the renormalized to the second order 0(e 2 ) effective interaction, eq. (|5"T|), as 



Af, 2 -Af' 2 X 
1/1 1 \ /. _,A!A 2 \ -I \ 2 ° 



2 Ui A 2 



K::\ = -e^ liA - hr + TT • l-e~W . e 



(52) 



l,\ 2 
MZ-M n 



where we have introduced 

9 



■ 





Ai _ A 1 


A - A2 - 


A 2 . 


PlP3 ~ 


~~ P+ (x'-x)P+ ' 


^P4P2 " P+ - (x> 


-x)P+i 




A = ^- ■ 

L - 1 P1-P2 P+ ) 


A - M 'o 





(53) 

(see Appendix B for the explicit definition of these quantities in the light-front frame). 

The expression eq. (|53| ) is written for the rescaled value of the potential V — > P +2 V, and 
the cutoff is defined in units of the total momentum P + , i.e. A — > p+ , with I = 1/A 2 . The spin 
structure of the interaction is carried by the matrix elements M 2 a, defined in Appendix B. 

We summarize the instantaneous interaction in both channels to the order 0(e 2 ), cf. 

fig- ©, 



_4e* 



-ir(ex) ^ C r r / 7*2 

*m.si :: 77+ _+ N9 SlS3 S2Si exp[ Pi. ^ 

(54) 



, ^ 4e 2 
T/;( e:c ) - _ X _ A _ exr>(-IA 2 ) 

v mst / + , +\ 2 u siS2 u s 3 S4 eA Pl iL -^PiPf) ' 

where we have used X^X^UXsiX^ = ^sis 3 ^s 2 s4 + ^sis 2 ^s 3 s 4 - For the rescaled potential in the 
light-front frame (fig. (§) and Appendix B eq. (|194|) and following) we thus have 



»W - (x _ x /)2 0-i-»°«-4 eX P 1 I A 2 



K msi,A ^ e 0sis 2 °s 3 S4 ex P ' \ A 2 / 



(55) 



where the notations of eqsgil4agil6 are implied. 

3.2.2 Renormalization issues 

As was discussed above after eq. (|TTD the commutator [t]^\ H eei ] also contributes to the self- 
energy term, giving rise to the renormalization of fermion and photon masses to the second 
order. The flow equation for the electron (photon) light-cone energy jT is 

j, ^ \J] i -^-ee-yl ^ self energy i V / 

where the matrix element is calculated between the dressed single electron (photon) states 
< p', s'\...\p, s >. We drop the finite part and define 5p\ = p~(l\)— < \H \ >. Integration over 
the finite range gives 

^ Q?S a (p)- (S£a(p)) 

n A P H 

that defines the cutoff dependent self energy 5Y>\(p). The mass correction and wave function 
renormalization constant correspondingly are given, cf. [El], as 



a - a - f h ^ r a) rr w AV (^a(p) -^ A (p)) 

0P\ ~0P A = < [T ',H en \ > se lfenergy dl = , (57) 



5m\ = p + Sp 



-<5£ A (m 2 ) 

(5* 



10 



The on-mass-shell condition is defined through the mass m in the free Hamiltonian H . 

We show further, that to the second order 0(e 2 ) the electron and photon masses and cor- 
responding wave function renormalization constants in the renormalized Hamiltonian vary in 
accordance with the result of 1-loop renormalization group equations. This can serve as evidence 
for the equivalence of the herediscussed method of flow equations and Wilson's renormaliza- 
tion scheme. Therefore we have rewritten the mass correction 5m\ through the self energy 
term, arising in 1-loop calculations of ordinary perturbative theory. The negative overall sign 
stems from our definition of the flow parameter, namely for Al > we are lowering the cutoff 
dl = -jzd\. 

We start with the bare cutoff mass m\ = m 2 + 5M A , where 8M A is the second order mass 
counterterm. According to eq. (^) the electron (photon) mass runs 

m 2 = m 2 _ [,5£ A ( m 2 ) _ 5S A (m 2 )] (59) 

defining, due to renormalizability, the counterterm 5M h = 5m\ = —5T l \(m 2 ) and the depen- 
dence of the renormalized mass on the cutoff A 

m\ = m 2 + 5m\ = m 2 — 5Y<\ . (60) 

We calculate explicitly the self-energy term. The electron energy correction contains several 
terms 

5p- x =< p\ s'\H - H \p, s >= (j2 SPxn) ■ S (3) (P - P%s> ■ (61) 

The first term comes from the commutator [r/(l), H eei ] 

SPi X = ~ < [T] {1) ,H eel ] >selfenergy dl' = ^ J (62) 

Jlx p 

it reads, cf. eq. ( |223| ) in Appendix D, 



2 r d 2 k ± dk+ 6(k+) .. + 1+ , 

xr( P - k, P , -k)r( P , P - k, k) l — — x (-r) . (63) 

p — k — [p — k) 

This term explicitly depends on the cutoff A through the /-function, that plays the role of a 
regulator in the loop integration 

^ = / p 2 M = exp|- 2 (^) 2 | . (64) 

Eq. ( |63D corresponds to the first diagram in fig. (0). 

Two instantaneous diagrams, the second and third in fig. (Q), contribute cutoff independent 
(constant) terms. Formally one can write 

drT - . dV in , st 

^- = < [v {2) , H ] > \ selfenergy =<00 + > -^- , (65) 

where < 00 + > stands for both the fermion and boson contraction (i.e. < b p b+ >= 0(p + ) and 
< a k a\ >= 8(k + ), respectively); and V™, st {l) = f pp >(l) V inst {l = 0) is defined in (). This gives 
rise to 

S PXn = < 66+ > -V^il) = < 66 + > -V inst (l = 0) (66) 

11 



for n = 2, 3. This means that 5p\ n defines together with Sp\i(l = 0) the initial condition for the 
total energy correction, eq. (|6l|) . 

In perturbative theory the instantaneous diagrams arise from normal-ordering Hamiltonian 
at Z = 0, and, in principle, must accompany the first diagram for any I. In what follows we use 
for the instantaneous terms the same regulator R, eq. (|6^), 



d 2 k ± dk + 0(fc H 



S P2X 

6P * X ~ ' 2(2vr)3 

We define the set of coordinates 



2(2vr) 3 k+ [p^ 



' 9(k^ 



k+] 
1 



x(-R) 



x 



k 



2\\p+ -k J - 



k+ 
p+ 

(xp + ,Xp ± + K ) 



\p" 



X 



-R) 



(67) 



(6* 



where p = {p + ,p ± ) is the external electron momentum. Then the electron self energy diagrams, 
fig. fl2|), cf. also eq. (|227|) in Appendix D, contribute 



p + 8p 



1A 



In 2 Jo 



dx / dn 2 , 



p- 



m 



«i + /(x) 



/(x) = xm 2 — x(l — x)p 2 



and 



2 
R 



2 + x - 



2m 2 



4 + /(x) 



+ 



m 



+ 



[1 



;r 



x(-i?) 



(69) 



p + 5p. 



2A 



^7T 2 JO 



dx / (iK 2 : 



„2 



x] [1 — x]J 



■R) 



p + Sp 



;U 



^7T 2 JO 



dx / dn i 



^7T 2 JO 



dx / dn 2 , 



2 



x 



(1 
-R): 



x 



X 



■fl) 



(70) 



for details we refer to Appendix D. Note, that the transformation in the integrals over x is 
performed before the regulator is taken into account 0. (In the second integral the electron 
momentum is replaced by the gluon one due to momentum conservation). The brackets '[ ]' 
denote the principle value prescription, defined later in eq. (|76|). 

The loop integral over k eqs. ( |69D and (|70|) contains two types of divergencies: UV in the 
transversal coordinate k, 1 - and IR in the longitudinal component k + . The physical value of 
mass must be IR-finite. We show, that the three relevant diagrams together in fact give an 
IR-finite value for the renormalized mass; this enables to determine counterterms independent 
of longitudinal momentum. In the wave function renormalization constant, however, the IR- 
singularity is still present. 

Define 

* = £ • ( 71 ) 
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where P = {P + , P ± ) is the positronium momentum, p the electron momentum. The transversal 
UV divergency is regularized through the unitary transformation done, i.e. by the regulator R, 
eq- © 

Rx = exp J - (^A 1 « 9(\% - |A p , fc |) , (72) 

where the cutoff is rescaled and defined in units of the positronium momentum P + , namely 
A — > a/2A 2 /P + , and A Ptk = p~ — k~ — (p — k)~ = K Pyk /p + . The rude approximation for the 
exponential through a ^-function changes the numerical coefficient within a few percent; never- 
theless it is useful to estimate the integrals in eqs. fl69|) and (|70D in this way analitically. From 
eq. (|7"2|) we have for the sum of intermediate (electron and photon) state momenta (the external 
electron is on- mass-shell p 2 = m 2 ) 

K L2 K ±2 + TTl 2 2 2 

TT + ~Ti i - ^ A <5i + m (73) 

[x\ [1 — x\ 

giving for the regulator 



R\ — 0(K Xmax K ) 0{K, XmaxJ 

K{l ax = x(l-x)X 2 5 1 -x 2 m 2 (74) 

and 0(K\max) leads to the additional condition for the longitudinal momentum 

^ _ ^ _ x max 

1 + m 2 l(\ 2 bi) 

implying that the singularity of the photon longitudinal momentum for x — > 1 is regularized 
by the function R\. This is the case due to the nonzero fermion mass present in eq. (|73|) for 
the intermediate state with (1 — x) longitudinal momentum. The IR-singularity when x — > 
is still present; it is treated by the principle value prescription || 

k+ = 2 \k+ + ieP+ + k+ -ieP+) ' ^ 

where e = 0+, and P + is the longitudinal part of the positronium momentum (used here as 
typical momentum in the problem being under discussion). This defines the bracket '| |' in 
eqs. (H) and (|70|) 

11/1 1 \ 

H " 2 [xT^ + V^ifJ ■ ( } 

Making use of both regularizations for transversal and longitudinal components, we have for 
the first diagram, eq. Q69|) , 



8m\ x = p + Sp 



IA _£ ( 3m ^ (**l£=!) + _^ (f^ + m .) _ 2A% ln («L^_ 

8 7r 2 I \ m 2 I X 2 di + m 2 \2 J \\ 2 b 1 + m 2 e 
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Note, that the third term has the mixing UV and IR divergencies. Combining the three relevant 
diagrams, fig. (|2]), and integrating with the common regulator, one obtains for the electron 
mass correction 



5m\ 
5m\ 



p + (Sp 1 + Sp 2 + 6p- 



3m In 



in z 



3/ \p 2 =m 2 



'\ 2 5 1 + m 2 



m z 



-cffi A (m 2 ) 
X 2 5im 2 
X 2 5i + m 2 



(79) 



The mass correction is IR-finite (that gives rise to IR-finite counterterms) and contains only a 
logarithmic UV-divergency. Namely, when \5\ —* A ^> m 



5m\ 



3e 2 2l A 2 
— -m In — - 



57T' 



77?/ 



(80) 



It is remarkable that we reproduce with the cutoff condition of eq. (|73|) the standard result of 
covariant perturbative theory calculations including its global factor 3/8. As was mentioned 
above, the difference in sign, as compared with the 1-loop renormalization group result, comes 
from scaling down from high to low energies in the method of flow equations. 

The similar regularization for the intermediate state momenta in the self-energy integrals, 
called global cutoff scheme, was introduced by W. M. Zhang and A. Harindranath ||. In our 
approach the UV-regularization, that defines the concrete form of the regulator R, arises nat- 
urally from the method of flow equations, namely from the unitary transformation performed, 
where the generator of the transformation is chosen as the commutator rj = [H Q , H]. Note also, 
that the regulator R, eq. (0), in general is independent of the electron momentum p + (rescaled 
cutoff X5\ — > A), and therefore is boost invariant. 

For the wave function renormalization constant, eq. ( J58f) , one has 



d5p 



dp~ 



in 2 Jo 



dn 



2 + x x(l — x)2m 2 



*l + ft) (4 + /(*)) 2 



x(-R) 



that together with the regulator R, eq. ([F^), results 
^ 2 ( X 2 5 x /3 n1 SA , 5i / 



1- 




m" 



2 In 



A 2 ^ 



1 — ' l i 

U X 2 5 1 +m 2 I 2 X 2 5 x + m 2 



In 



1 



Si 



X 2 5, 



2 A 2 ^ + m 2 
As XSi — >• A ^> m the function F tends to a constant 



In — I + F In 



-2 + 2 



X 2 5 x 



X 2 5 1 



X 2 5\ + m 2 X 2 Si + m 2 
lnx 



dx- 



x — 1 



(82) 



|A>m 



c 



3 

2 



Therefore, by dropping the finite part, we obtain 

.2 r A 2 /3 



l 



57T' 



In 



m' 



2 In 



7T" 

y 



+ ln 




(83) 



(84) 



where we have rescaled ■# — > e. The electron wave function renormalization constant contains 
logarithmic UV and IR divergencies mixed, together with pure logarothmic IR divergencies. We 
mention, that the value of Z 2 is not sensitive to what kind of regulator is applied; the same 
result for Z 2 was obtained with another choice of regulator || . 
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We proceed with renormalization to the second order in the photon sector. The diagrams 
that contribute to the photon self energy are shown in fig. (||]). The commutator fr/ 1 -*, H en ] , 
corresponding to the first diagram, gives rise to (cf. eq. ( |235| ) in Appendix D) 

^ - bV7w s(fcW - fc+) (85) 



xTr 
where momenta are given in fig. (Rt) , and the regulator is 



T\k, k - q, q)T j (k - q, k, -q)\ ^_ _ fc _ _ _ X (-R) 



A 



q.k 



2" 



*A = &,A = exp|-2(=3£) | . (86) 

In full analogy with the electron self energy this also defines the regulator for the second diagram 
with the instantaneous interaction, see fig. (^), 

°* s " = wf I w" (fc+) Tr(ffV) (f^T - b^Ff) x ( - fl) ' (87) 

We define the set of coordinates 

(q-k)+ 



= x 
q^ 

k = ((l- X )g+ ) (l_ X )g ± + K ± ) 

(q-k) = (xq + , xq 1 - K L ) , (88) 

where q = {q + ,q ± ) is the external photon momentum. Then two diagrams contribute (for 
details see Appendix D, eq. (|239|) ): 



q + Sq 1 = —r—z dx d,K 2 ± 
Stt^ Jo J 



x { 4T7M ^ ~ 2x + + 4T7W + (" 2 + mpM } x { ~ R) 



f(x) = m — x(l — x)q 

e 2 r°° , r , , / I I 



3tT 2 



dx J dKl {iT-z\-—x) x{ - R) 

f dx /d«lA x (-R) . (89) 

Jo J \x\ 



Note, that the transformation in the second integral is done before the regularization (by reg- 
ulator the R) is performed ||. 

Making use of the same approximation for the regulator as in the electron sector, we obtain 
for the sum of intermediate (two electron) state momenta 

k 2 ± + m 2 k 2 ± + m 2 2 

1 <Ad 2 

x 1 — x 

^ = J^ , (90) 
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where the photon is put on mass-shell q 2 = and the rescaled cutoff A 
used. This condition means for the transversal integration 



K 



Rx 

J-2 
Xmax 



6{k 



±2 
Xmax 



— K~ 



0(k 



Xmax J 



x(l — x)\ 2 S2 — m 2 



v / 2A 2 /P+ has been 



(91) 



and for the longitudinal integration 



a?i < x < %2 




1 — r 


2 
771 


2 ' 


A 2 <5 2 


1 + r 

%2 = „ « 


2 



'1 



4m 2 
A 2 ^ 



f92) 



where the approximate value is when m< A. This shows that the condition of eq. fl9U| ) for two 
electrons with masses m removes the light-front infrared singularities from x — > and x — > 1. 
Thus, both UV and IR divergencies are regularized by the regulator R, eq. (|9T|). 
The mass correction arising from the first diagram, eq. (B9l), is 



e 2 2 



1A 8tt 2 3 



= A 2 o^ 



4m 
A 2 ^ 



2\3/2 



(93) 



Combining together both diagrams with the same regulator, eq. (p9|), we obtain 



<5m 2 



57T 



' \2jt 8 2 o 2i 1+r 

-A 02 r m r — 2m In 

1 3 1-r 



(94) 



where r is defined in eq. (^). The result shows that the mass correction involves the quadratic 
and logarithmic UV divergencies, i.e. as A5 2 —* A ~^> m 



5m\ 



in 



-J S -A> 

2 I 3 



A 2 " 
2m 2 In — - 
m 2 , 



(95) 



The wave function renormalization constant is defined through 



d5q 



dq~ 



in 2 Jo 



1 r ( 2x 2 

dx / dn 2 , 



q 2 =0 >-"' "u - ^ k ± + j(x) 

that, with the regulator R, eq. (0), results 

,2 / 2 ^ 1 + r 



2x + 1 2m 2 x(l — x) 



mr- 



57T" 



, . ,10 8 wr 

In 1 r -\ 

3 1-r 9 9\ 2 5 2 



x 



-R) 



(96) 



g 2 =0 



(97) 



The photon wave function renormalization constant contains only logarithmic UV divergency, 
indeed as Ao~ 2 —> A ;§> m 



1 



in z 



2 A 2 
v 3 m 2 ^ 



(9* 



and is free of IR divergencies (as is expected from the form of the regulator R, eq. (|90] 
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3.3 Renormalized theory to the order e 2 (at the scale A) and Feyn- 
mann rules 

We have completed the renormalization of light front QED (LFQED) to the second order, what 
enables us to formulate the diagrammatic rules for the perturbative expansions in cq. We formu- 
late the Feynman rules (in light front frame) for the renormalized x + -ordered LF Hamiltonian. 
All matrix elements of the renormalized Hamiltonian H r , namely the interacting part of H r 
with the corresponding LFQED vertices, are listed in fig. ([!]) together with their diagrams. The 
diagrammatic rules are obtained by direct calculation of matrix elements between free particle 
states. 

The two last diagrams in fig. ([!]) correspond to the generated interactions, arising from 
the renormalized Hamiltonian considered in different sectors. Note, that the existence of the 
|ee7 >- vertex in the renormalized Hamiltonian at any finite cutoff A prevents the generated 
interaction to have small energy denominators 4-, 4- (the |ee7 >-vertex appears through the 
exponential factor in the generated interaction). We mention, that the problem of small energy 
denominators in generated terms also is solved in the similarity scheme of Glazek and Wilson 

B- ' , 

The two-component LF theory, introduced by Zhang and Harindranath ||, as compared 
with the four- component formalism of Lepage and Brodsky is formulated purely in terms of 
physical degrees of freedom, so that each term corresponds to a real dynamical process. This 
means, for instance, that the first of the two generated term diagrams must be taken into 
account in low energy ee-scattering, and the second one describes Compton scattering at low 
energies. Further, we use the instantaneous and generated interactions, namely the second and 
the forth diagrams in ee sector, together with the perturbative theory contribution to the order 
0(e 2 ) of eej vertices (in the same sector), to calculate the mass of the positronium bound state. 
The rules to write the expression of perturbative expansions from diagrams are given in ||. 

Note, that the similarity function f PiPj \ in each term restricts the energy differences between 
initial and final states to be below the UV cutoff A, and thus restricts the renormalized Hamil- 
tonian to act in the low energy sector. Therefore, for the UV cutoff A << A, we associate the 
renormalized Hamiltonian with the effective Hamiltonian descibing the physics of low energy. 

Another aspect of the flow equation method is connected with renormalization group. 
Namely, the function f PiPf \ in the ee7-vertex plays the role of UV (and partialy IR) regulator 
in the self energy integrals (see before), so that the regularization prescription of divergent inte- 
grals originates from the method of flow equations itself. Moreover, as we will further show, the 
energy correction (i.e. mass correction and wave function renormalization constant), obtained 
from the flow equation method, coincide up to the overall sign with the 1-loop renormalization 
group result. This is the remarkable result, indicating the equivalence of flow equations and 
Wilson's renormalization. 

At last we mention that in the diagrammatic rules given in fig. (|I|) we explicitly write the 
dependence of the electron (photon) mass on the cutoff, namely 

m\ = ml- 5S A , (99) 

where 5T,\ is the self energy term to the order 0(e 2 ), (m = for a photon). Whereas we drop 
the subscript A for the polarization vectors e and spinors x- To the next order 0(e 3 ) one has 
e x = elil + 0{e\)). 

In the next section we use these diagrammatic rules of LFQED renormalized Hamiltonian, 
fig. (p, to calculate perturbatively electron-positron interaction, namely the contribution from 
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perturbative photon exchange and electron (photon) self energy correction, which latter enables 
to find the corresponding physical masses. 

4 Positronium's fine structure 

4.1 LF Perturbative theory 

The scattering |ee > states are also needed in bound state calculations. Using the propagator 
techniques we include these states where required. We exploit the perturbative theory in the 
coupling constant e, using the Feynman rules of the renormalized theory fig. ([!]). 

The first order renormalized ee'j- vertex f PiPf H^ contributes to the second order to the 
\ee > interaction term and to the electron (photon) mass renormalization. Physically, it is the 
perturbative photon exchange (photon emission and absorbtion in the case of electron mass 
renormalization), with the energy widths of the photon restricted by the function f PiPf ,\- 

4.1.1 The electron-positron interaction 

According to the light-front Feynman rules the perturbative photon exchange gives rise to the 
following second order |ee > interaction in the exchange channel 



V(l) = g x {l)g 2 {l)M 2ii ■ I ""> - l - + 7 q J - l - 1 , (100) 

I <r Pi - p k {-q + > Pi -Pk ) 

where Qi stands schematically for the coupling constants in both vertices, namely g plP2 x = 
ef plP2 x and M 2 u defines the spin structure of the interaction, coming from the corresponding 
structure of the eej- vertex; and the two terms in the curly brackets represent two different x + 
(time) orderings of the photon exchange with the momentum q, giving rise to the two different 
intermediate states with momenta p^, Pi corresponds to the initial state. Explicitly one has 
In the exchange channel 

vP(i) = - e 2( / _ p4i _ p2 (/) x +r(-p 4 ,-p 2 ,-g) Xs - 4 )(/ piiP3 (/) x +r(p 1 ,p 3 ,g)x, 1 ) (101) 

'9(p+-pt) 1 , 0(p+-pt) 1 



x 



°q,pi-P3 



(Pi-Pt) Pi -Ps -Pi -9 ipt-Pi) Pt -Pi -Pa +1 

with the initial state momentum pi = P = (P + , P ) and momentum transfer q = p\ — ps, and 



P ±2 + Ml 
~P+ 



P- = ^— -^ , (102) 



where Mjy is the mass of positronium bound state. In the light-front frame 

- ( p ~ -Pa -P2 ~ {Pi-Pz)~)°{Pi -Pi) = (P~ -Pi -Pa + (pi - Ps)~)9(pt - Pi 

A 3 



P+(x - x'] 
holds, giving for the rescaled potential V —> P +2 V rise to 

r(ex) ___ 7 at 1 m ^( (Af + A|j 



(103) 



V±Z = -JKj-exp [- X i ) . ( 104 ) 
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where iVj is defined in eq. (0) and 



A 3 = (k ± - k' ± ) 2 + ~(z - x')A + \x- x'\ Q(M 2 + M >2 ) - M 2 N 



1 1\ ,., 2 , 2 , /l 



(105) 

A = (** fm 2 ) (-_--) + (*? + m a ) f ^ - -J-^) , 

VI — a; 2/ Var I-17 

Here the cutoff A is defined in units of P + . 

Because of the absence of Z-graphen in light-front formalism (corresponding to negative 
p + ), only one term contribute to the annihilation channel, namely 

Vfe°(0 = e 2 {f pl ^ P2 {l)xf 2 T l {p 1 , -p 2 , q)Xs 1 ){f-p i , P3 {l)xf s ^{-P^P3, -q)Xs 4 ) 

$q, P l+Pl > ( 106 ) 



(pi + pi) Pi -or 



where p i = P and the momentum transfer is q = p\ + p 2 - This gives rise for the rescaled 
potential V — * P +2 V in the light-front frame, to the expression 



1 j (jg + M '*' 

" 2 mT xp V 



V& = A^«p-^W, don 



£ iV 



where ./V2 and the variables Ai, A2 and Mq , M' are defined in eq. 

4.1.2 Mass renormalization 

Following light-cone rules the perturbative energy correction of the electron with momentum 
p, coming from the emission and absorption of a photon with momentum k, is 

6 ^ x = J 2(27r) 3 k^~ 6<KP+ ~ k+ ) g P~ k ^ r ^' p ~ k > p > - k )9 P , P -k,\T\(p,p- k, k) 

V-^-(p-*>- • (108) 

where g ee7 -coupling constant restricts the energy of the photon. Making use of the explicit form 
for the coupling, one has 

«a - e2 /w^ (p+ - fc+) (109) 

xrl(p - k,p, -k) T\(p,p - k, fc) — p-L— — - x (R) , 

where R = f 2 kX plays the role of regulator. This expression coincide up to the overall sign with 
the energy correction obtained in the previous section from the flow equations method. 

Two instantaneous diagramms, arising from the normal-ordering Hamiltonian, must be 
added to the first term with the same regulator R. Then the full perturbative energy correction 
$P\ = SPlx + &Pi\ + SPsx is 

% = -fa ( 11Q ) 
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where 5p x is defined in eq. fl6Tp. This means for the perturbative mass correction 



5m p x T2 = 5S A (HI) 



and the self-energy term £E A is given in eq. (|79|) . 

We combine the renormalized to the second order mass, eq. (pOD, and the perturbative 
correction, eq. (|111|) , to obtain the total physical mass to the order 0(e 2 ) 



m 2 e = m{ + 5m 2 = (m 2 + <5£ A ) - <5£ A = wf + 0(e 4 ) . (112) 

This means, that to the second order 0(e 2 ) the physical electron mass is, up to a finite part, 
equal to the bare electron mass, that stands in the free (canonical) Hamiltonian. 
Along the same line one can do for the photon mass. 

4.2 Bound state perturbative theory 

In this subsection we define bound state perturbative theory (BSPT). 

First introduce instead of the front parametrization, used before for the single-particle 
momenta, fig. (^), the instant form 

p Ul = (xP + ,xP^ + k L ,Pi) — ^ Pul = (k z ,k ± ,p1) = (pi, Ex) 
p 2{1 = ((1 - x)P+, (1 - x)P ± - k^ ,p 2 ) — > p 2fl = (-k z , -k L ,p Q 2 ) = (p 2 , E 2 ) 

E { = y/p 2 + m 2 , i = l,2 



x 



Ex+k z _ 
E!+E 2 2 



h i + 



yp 2 + r 



113) 



and for the momenta P3,P4 the same, but with prime over x,k z ,k x ; here x is the light-front 
fraction of the electron momentum, and J(p) is the Jacobian of the transformation: 

. x Q/OC fv I — r" Til . , 

[P) = dk~ z = 2(j^ + m2)3/2 • ^ 114 ^ 



The new coordinate system corresponds to the center of mass (cm.) frame, i.e. P = 0. Note, 
that in both sytems electron and positron are moving along paths on the light-cone (against 
each other in the instant case P3 = k z \ + k z2 = 0), while the connection between the two frames 
is obtained by boost. 

The light-front bound state equation (see later) is boost and frame invariant; therefore it 
can be solved in the cm. frame, producing the same spectrum as in the front form. We 
shall see that the electron-positron potential, arising from the renormalized to the second order 
Hamiltonian, is drastically simplified in the cm. instant frame. Also, this system is usefull 
in practical aspect: in making obvious the rotational symmetry, restored in the nonrelativistic 
limit and manifest in the spectrum. 

We stress an important consequence of this change of coordinates. Since the total light- 
front momentum P + = Pq + P3 is conserved (i.e. both 'in' \x, k x > and 'out' \x', k ,J - > Fock 
states have the same P + ) and P3 = in instant form, we have conservation of the total energy 
Po = Pq, or, in other words, E\ + E 2 = £3 + £4. For the 3 — d momentum this means 

f = p' 2 . (115) 
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This condition we will further exploite in what follows. 

Define now BSPT. We choose the leading order electron-positron potential in such a form 
to simplify positronium bound state calculations. This means, that this potential contributes 
the leading order term to the positronium mass, and perturbative theory with respect to the 
difference between the total second order |ee > interaction, calculated before with the renor- 
malized Hamiltonian H r to 0(e 2 ), and the leading order potential converges. This scheme we 
call BSPT. Surely, our choice is motivated by the form of the renormalized to the second order 
interaction to insure convergence of BSPT. 

We define the second order renormalized electron-positron potential (e(3)e(4)| V cou i |e(l)e(2)) 
to the leading order of BSPT in the form of pure perturbative one photon exchange, explicitly 
as the Coulomb interaction 

16e 2 m 2 16e 2 m 2 

coul = ~(k±- WJ* + (k z - W = ~Jf^Y ' ( } 

This means that the corresponding leading order Hamilton operator in the \ee > sector is 

H^ = h + V coul , (117) 

where h is the free part, defined in eq. (0). The wave functions are given as the solution of 
Schrodinger equation 

H^\i/j n (P) >= E n \^Pn(P) > , (118) 

where P is the positronium momentum, and the eigenvalues and eigenfunctions for the positro- 
nium bound state are defined in standard way of light front frame 



E 



N 



P+ 



IMP) >=E S1S2 I P1P2 ^pipi2(27cf5^(P-p 1 p 2 )$ N (xk ± s 1 s 2 )bf i (p 1 )df 2 (p 2 )\0> 

f d?k f 1 dx ~ ~ 

^siS2 2(27r)3 $ ^(^±SiS 2 ) & N (xk±S 1 S 2 ) = 5 N N' (1,19) 

Mn stands for the leading order mass of positronium. Combining the definitions for the wave 
function and the energy with the Schrodinger equation, we obtain 



M\ 



k'l + m 2 1 ~ , , t-^ J d 2 k±J^ dx 



2 _ 

v x'(l - x') 



$ N (x'k' ± S3S 4 ) = J2 o/o x° 3 V coul ®N (xk±SiS 2 ) , (120) 

sis 2 2,\2m) 



or, after change of coordinates according to eq. (|113| ) 



_ f d 3 p\ J(p)J(p') 
M 2 N - 4GT + m 2 ) ) $ N tfs 3 St) = V ■ 1 \ 3 V coul (p,p>) <5> N (ps lS2 ) , (121] 



S1S2 

where the wave function was redefined to have the norm 



J2 I d 3 p$* N (p Sl s 2 ) & N (p Sl s 2 ) = <W . (122) 

S1S2 J 



We aim to obtain the nonrelativistic Schrodinger equation for positronium. Note, that in the 



rr~ 



nonrelativistic limit ^ « 1 we have 



'mm M i (i - t»&^i 

M N = {2m + B N f « 4m 2 + 4mB$ , (123) 
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where we have introduced the leading order binding energy B N . Then to the leading order the 
bound state equation for positronium is 

% - B^ $ N tfs 3 St) = -E S1S2 

Making use of the explicit form for the Coulomb potential, eq. ( |116| ), we obtain the equation 
that determines the leading order bound state wave function: 

*(p') = ^L f dp ^ $ (p) (125) 

^ ' 2tt 2 J (p-p>) 2 M 

with 

$tv = $Hs.,M,(fis 3 S4) = $„(??) S SeS3 S ssst . (126) 

This is the standard nonrelativistic Schrodinger equation for positronium. Its solution is 
characterized by fi = (n,l,m), the usual principal and angular momentum quantum numbers. 
The wave functions are given through the hyperspherical harmonics 

(e 2 + V 2 ) 2 

Yp = Y n ,l,m = fn,l(u)Yl,m(6,<P) 

ma 2 ma 

B " = '!*■ e " = *7 < 127 > 

and for the binding energy one has the standard nonrelativistic expression for positronium 
bound state to 0(e 2 ). For sake of completeness we write the coordinates used in the solution 

(e 2 n = -mB N ,p) — > (u ,u) 
i -*! 

U = COSCJ = 

ei+P 

u = prsmw = , (128) 

IP] ei + p 2 

but, for details, refer to |J. 

The electron-positron interaction arising from the renormalized to the 0(e 2 ) Hamiltonian is 
given as a sum of two contributions from exchange and annihilation channels V = V exc h + V ann 
(see explicitly later). We introduce the potential, arising in the nonrelativistic Schrodinger 
equation, eq. ( |125| ), 



<J(p)J(p>) l 
V r (ps 3 S4;psis 2 ) = Jim — — — (V exch + V ann ) . (129) 

Then we define BSPT with respect to the difference 

~ - a 1 

SV = V[p's 3 S4;psis 2 ) - ("2^2) /~_ -^ 'm 5 **** ' ( 130 ) 

where the leading order contribution is defined in eq. ( |127| ). Note, that, to define the Coulomb 
potential, we have taken only the first term of the Jacobian's J(p) nonrelativistic expansion, 
i.e. the ee interaction to the leading order of BSPT. 
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In what follows we use the matrix elements of 8V, defined as 



< ®nim\5V\<5> n i m >= J d 6 pd 6 p®* nlm (p)8V® n i m {y! 
where $ n ; m are the Coulomb wave functions given above. 



(131) 



4.3 The renormalized electron-positron interaction in light-front and 
instant form frames 

We summarize all together the second order |ee> interaction in the exchange and annihilation 
channels, i.e., written in the front form frame: 



V, 



exch — V\ + V — V\ + V\ + V\ 



Vann = V^ + V™ = ^ »"* + ^ ^ + V{ PT , (132) 

where the generated, instantaneous and perturbative theory interactions (rescaled, i.e. V — > 
P +2 V) are given correspondingly in the exchange channel 
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in the annihilation channel 
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where in the light-front frame, see eq. (|194f ) in Appendix B, 



-<'1,A — "SlS3"S2S4-' 1 ■-'2 ' "S1S2"S1S3"S2S4^^^ 



(x — x') 



A2 



arx'(l — x)(l — a;') 



+zmv2(x / — x) 



r _ r _^_7^-L -J- _|_ X X f^ T" 1 - 



S 2 



iV. 



2. A — SlS2 S 3 S4 J 3 • i 4 + S iS2"s 3 S4"siS3^ m ^/Q _ .Ui _ „A 

1 ■ /ox X Sl T*-L r -L _ A A g 3 ^_|_ _|_* 

'/"\ - s3f4 sis2 x d _ x ) 3 s i 3S4 i52 x'(i — x'j 4 4 



V2 



s,z 



(135) 
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and 
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T' 
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with the definitions 



Ax 



(xk' ± — x'k±) 2 + m 2 (x — x') 



[1-x) 

A2 
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x'(l — x') 
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P x ) 2 + M 2 
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1 — x x 
v P = (P+ P 



1 



x' 1 — X 1 , 

M N = 2m + B N 



(136) 



Note, that the rescaled potential, eq. ( |132| ), does not depend on the total momentum P + , i.e. 
is invariant under light-front boosts. 

The generated interaction plays an important role, namely it insures the absence of collinear 
divergencies in the renormalized interaction, which are associated with the limit x — > x'. This 
is true for any cutoff A. In fact, in the limit when x tends to x' one has 



Ai ~ A 2 — ► (k ± - k ± ) 2 
A A ^ ~ k '^ 

Ai ~ A 2 > ; 

X — X 
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This gives rise to 

V* xch {x -> x') 



e 2 N, x 



(k± - k' ± 



_ r gen 
2 ex 



4e 2 



[x — x 



c inst S S 

l\2 ex U S!S3 U S2S4, 



( <-< )' 



'1381 



where, in this limit, the perturbative part vanishes. Making use of eq. ( |135| ) we obtain for the 
divergent part DV\- ■ ■] 



S2S4 



VV[N l:X (x^x')}=T 1 ± T 2 ± S SlS3 S. 

VV[T{{x - x')} = -2 %— k -4- ; VV[T 2 {x - x 1 )} = 2^ ~ k±)% 



[x — X' 



(x — X 



(139) 



This results in the exact cancellation of the collinear divergencies present initially in the in- 
stantaneous interaction by the generated term: 



VV[V^ xch (x -► x')} = 



(140) 
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Note, that no approximation for the ee interaction was made untill now. 

We rewrite both exchange and annihilation channel contributions in the instant frame 



V — V exc h + V Q 



exch ~T~ V ann 
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+ (V 8 si s 3 5 S 
where quantities in the instant form are defined as follows 



1411 



x 



Ml 




2 y \/p 2 + rri z j 



M A =4(f + m 2 



(142) 



for the other quantities, defined in eqs. ( |135| ) and ( |136| ) the substitution x(k z ),x'(k z ) is to be 
done. The symbols cf^ n and so on were introduced to indicate the origin of the different terms 
(here generated interaction coming from the exchange channel), all c = 1. We stress once 
more, that at the moment no approximation was done, therefore the expression of eq. ( |141| ) 
can be considered as the exact electron-positron interaction, arising from the renormalized to 
the second order Hamiltonian, that produces to the leading order the positronium mass M^ 
and the binding energy B^. In fact, it will be shown, that the second order potential gives rise 
also to the correct mass splitting, that corresponds to the next to leading order contribution. 

The expression of eq. (|141|) has a transparant form. First notice, that Ai and A2 descibe the 
energy differences in two corresponding ee7-vertices appearing in the ee interaction. Then the 
generated interaction (cf^ n ) contributes mainly hard photon exchanges ->4 ~ -rf >> 1, while 
the term arising from perturbative theory c^J gives rise to soft photon exchanges. Though 
the renormalized interaction generally descibes low energy physics, namely the renormalized 
Hamiltonian acts in the space where the energy differences are restricted by A, the information 
on the high energy sector is accumulated in the generated interaction, making possible to 
interpolate between two sectors. This means, that the sum of both terms in eq. (|141|) recovers 
the whole range of photon energies. The same is true for the annihilation channel. 

Two limits are of interest. The trivial limit A — > A reproduces the bare interaction, used 
as initial condition for the renormalized interaction as A tends to infinity. In the opposite limit 



A «. m 



(143) 



nonrelativistic physics is affected (see later). The limit A — > can be performed explicitly, 
corresponding to the complete elimination of the ee7-vertex present in the initial Hamiltonian. 
When A tends to zero the ee interaction is governed by generated and instantaneous terms. 
The fact that in this limit the potential is well defined is the direct consequence of the instant 
form frame. 

The expression of eq. ( |141| ) for the renormalized to the second order interaction, written 
in the instant form frame, can be used as a suitable form to obtain the positronium spectrum 
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numerically. We further proceed in another direction, namely we make the nonrelativistic 
approximation (NR) for the electron (positron) momentum 



\P\ 



m 



0(a) < 1 , 



that enables to perform calculations for positronium mass splitting analitically. 
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4.4 The nonrelativistic approximation 



In the nonrelativistic approximation we obtain for the ee interaction, eq. (|141 
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where the energy denominators and exponential factors were simplified using 



x — x 



^2 K 

2m 



1 + 



-o 

P 



2m 2 



+ 0[m 2 



P 



2 P 



m 



A 1 = A 2 = (p-piy + 0[m 2 [^- 



711 



\2 , I™ '1/^2 B/f 2\ , /-> / _2 / P 



A 3 = (p - p') 2 + |x - x \(M 2 -M 2 N ) + 0\m 2 



in 



Ai = A 2 



2m(p' — p)' 



1+0|[£ 
m 



A 



2m(p' — pf 

(K - k z ) 



Ml = Am 2 + 0[m 21 '' 

m 



M% = Am 2 + 47715^ + O ( m 



5^ 2N 



2 / -pjv 

in 



Am 2 + 4m5 



(o) 

JV j 



and the explicit expression of Jacobian for the coordinate change is 



J(p)W = ^ 



!_Jl + (4,^ 
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(146) 
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having introduced the leading order binding energy B N . Making use of its nonrelativistic 
approximation B N /m « 1 we have for the interaction 
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where (c 9en , c PT ) showes that both term, generated and perturbative interactions, contribute 
to corresponding term (we remember that all c = 1). 

The remarkable feature of the part of interaction standing in the first bracket is that it 
does not depend on the UV cutoff A. The next term in the second bracket arises from the 
perturbative photon exchange and has the typical 'energy shell' structure for the relativistic 
effects, namely these terms are important when A >> m. Further we calculate the ground state 
positronium mass and therefore restrict the cutoff to be in the nonrelativistic domain 



A «. m 



(149) 



where the second term in eq. ( |148[ ) vanishes and we are left with the following form for the 
renormalized ee interaction in the nonrelativistic approximation: 
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We perform the nonrelativistic expansion of the factors Ni and N 2 appearing in the interaction. 
The term N\ contributes in V to the order 
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Whereas the term A^ contributes to V to the order 
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In these formulas we have used [k ± , fcj_ 
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following variables have been introduced 

q± = k\-k±. 



?ij 3 and e\ = — -k=(s,i); also the 
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We leave aside for the future work the analysis of the expressions for iVi and A^, where also in 
this form some terms can be identified as spin-orbit and spin-spin interactions in the transverse 
plane and in longitudinal (z) direction. 

Instead we follow |J, where an analogous calculation of singlet-triplet ground state mass 
splitting of positronium was performed in the similarity scheme. This means, that we can, 
except for the leading order term 0(1), drop in Ni the part diagonal in spin space. Also the 
terms of the type / = k± y k z , k x ± y k z , k\k\ do not contribute to the ground state mass 
splitting, since 
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averaging over directions, gives zero. 

We obtain for the ee-potential to the leading order 0(1) of nonrelativistic expansion 
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Remembering q = p'—p, Fourier transformation to the coordinate space with respect to q has 
been performed in the last expression. To the leading order of NR expansion we have reproduced 
the Coulomb potential, defined before as the leading order of BSPT. Note, this is true for any 
UV cutoff within the nonrelativistic range A « m. 

We combine this expression with the kinetic term from the Schrodinger equation, eq. ( |125|) , 
and write it in the form 

- f 1 + ^P) f + V(r) . (155) 

m \ 2m I 
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Here the potential V(r) plays a different role in the two terms. In the first term, corresponding 
to kinetic energy, it generates an effective mass of the electron, which depends on the relative 
position and manifests the non- locality of the interaction. The second term is the usual potential 
energy, in our case, the Coulomb interaction. 

The energy of the Coulomb level with quantum numbers (nlm) is standard 

Ml =< $ nlm \V^ nlm >= J d 3 pd 3 p'^ nlm (p) V® $ nlm (P) = -~ , (156) 

where the Coulomb wave functions $ n z m were defined in eq. ( |127|) . We have used in eq. ( |156j ) 
the following representation 
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and also orthogonality of the hyperspherical harmonics 

^<mY;Yj = 8^. (158) 

More details can be found in |J. 
The next to leading order O f %- ) 

2(27r) 3 4m2mV q 2 ) 

x (8m(iql - s 1 q y ± )5 SlS3 5 S2S4 - 8m(iql - s 2 q y ± )S SlS3 5 S2 s 4 ) (159) 

contributes (because of the spin structure) to the second order of BSPT: 
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2 ~h Ml -Ml ■ {m 

The order 0((^) 2 )(cf. remark after eq. «) is 
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and contributes to the first order of BSPT: 

5M 1 2 =<$ 10 o|^ 2 )|$ 10 o> • (162) 

Both contributions were calculated in |J with the result 

5M 2 = 5M 2 + 5M 2 

< 1\5M 2 \1 > = -—ma 4 
1 ' 12 

<2\5M 2 \2> = <3\5M 2 \3>=<4\8M 2 \4>=-ma 4 , (163) 
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where the eigenvectors in spin space are defined as follows: 

|2 >= -7= (|+ -> +|- +>) , |3 >=|-->, |4 >=|++>. (164) 

v2 

Making use of the relation between Coulomb energy units and Ryd = |ma 2 we have the 
standard result for the singlet-triplet mass splitting for positronium, |a 2 Ryd. The degeneracy 
of the triplet ground state n = 1 reflects the rotational invariance, manifest in the system in 
nonrelativistic approximation. 

5 Conclusion 

Basing on the similarity scheme together with the flow equation method we derive the renormal- 
ized to the second order 0(e 2 ) LFQED Hamiltonian. It has, in its turn, two aspects considered 
further in the work. 

The first point is related with renormalization group. It turnes out, that electron (photon) 
mass correction and wave function renormalization constant vary with UV cutoff in accordance 
to 1-loop renormalization group equations. This indicates an intimate connection between 
Wilson's renormalization and the flow equation method. The same equivalence to Wilson's 
renormalization was shown for the projection operator techniques (Bloch-Feshbach formalism) 

0- 

In the case of LFQED two types of divergencies arise: UV divergencies, associated with 
a large transverse momentum and severe IR divergencies, originating from light-front gauge 
singularities and present in light-front Hamiltonian calculations for QED. Flow equation method 
deliver in divergent integrals straightforwardly the regulator, that regularize UV and partially 
IR divergencies. It was shown in the work, that gauge invariant calculations of divergent terms 
give rise to IR finite results. Explicitly, the IR divergencies were removed from the electron 
(photon) energy correction, when all diagramms to the second order in the corresponding sectors 
(arising from flow equations and normal-ordering Hamiltonian) were considered. This enables 
to choose IR-independent mass counterterms and insures IR-finite physical masses. 

To complete the renormalization group analysis the third order (coupling constant renor- 
malization) and Ward identities must be considered in the techniques discussed. 

The second aspect is connected with the low-energy sector of the theory and bound state 
calculations. The renormalized Hamiltonian acts in the space where the energy differences be- 
tween initial and final states are restricted by the UV cutoff, and therefore it can be interpreted 
as an effective Hamiltonian descibing the low-energy physics. 

Making use of diagrammatic rules for the renormalized LFQED Hamiltonian, we obtain 
the electron-positron interaction to the second order. The generated interaction, arising from 
elimination of ee7-vertex, insures at any cutoff A the absence of collinear divergencies in the 
ee-interaction. 

The ee-interaction possesses the explicit and implicit, through the running mass and cou- 
pling, cutoff dependence. We aimed in the work to get rid of explicit A-dependence, performing 
the limit A — > or making use of the nonrelativistic approximation, to obtain the physical 
results for the spectrum. The instant form frame was used, that manifest several advantages as 
compared with the light-front one. First, the complete elimination of the ee7~ vertex is possible, 
that corresponds to the limit A — > where the ee-interaction obviously does not depend any 
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more on the cutoff A. In this limit the interaction is governed by generated and instantaneous 
terms. Note, the light-front bound state equation is boost and frame invariant, giving rise to 
the same spectrum in both frames. 

To calculate the positronium spectrum analyticaly the nonrelativistic approximation with 
respect to the electron momentum — = 0(a) and binding energy — - = 0(a) was performed. 
The instant form enables then to obtain the cutoff independent ee-interaction as long as the 
UV cutoff is restricted to be in the nonrelativistic domain, namely A « m. The leading 
order (in nonrelativistic expansion) of this interaction describes 3d Coulomb interaction, giving 
rise to standard Coulomb energy levels and wave functions. Note, that in the nonrelativistic 
approximation it occurs exact cancellation of instantaneous interaction by the leading order 
generated and perturbative theory terms, giving rise to the 3d-Coulomb interaction. The next 
order terms define the spin structure of the ee-interaction. The same result for the spin- 
dependent terms was obtained in |J, where the ground state singlet-triplet mass splitting for 
positronium was calculated to be ^a 2 Ryd. The mass degeneracy of the triplet ground state 
manifests the rotational invariance, maintained at the order a 4 . 

The correct results obtained for the positronium mass spectrum in nonrelalivistic approxima- 
tion are encouraging. This makes possible to consider an exact expression for the ee-interaction, 
written in instant form frame, as a basis for future numerical calculations. 
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A Renormalization scheme of similarity unitary trans- 
formation 

Here we discuss briefly the renormalization scheme presented by Glazek and Wilson ||, called 
similarity unitary transformation, and compare it with the flow equations by Wegner. 

Both approaches aim to bring the Hamiltonian to the most diagonal form, explicitly only 
the matrix elements between the free states with 

|Ay|<A (165) 

and Aij — Ei — Ej, are present in the renormalized Hamiltonian. For this purposes the con- 
tinuous unitary transformation must be performed to preserve unchanged the spectrum (eigen- 
values) of the initial bare cutoff Hamiltonian. The demand of diagonal structure does not 
define completly the generator of the transformation. This freedom is used to eliminate small 
energy denominators in the final renormalized Hamiltonian. This results in a system of two 
self-consistent non-linear differential equations for the Hamiltonian H(l) and the generator of 
the transformation T](l). The dependence on the continuous flow parameter I in the flow equa- 
tions by Wegner is replaced by the cutoff dependence A in the similarity approach, with the 
connection 

I = 1/A 2 (166) 

in the renormalized Hamiltonians. 

The difference of the two methods manifests the residual freedom in the choice of the direc- 
tion of the infinitesimal rotation, actually defining how fast the non-diagonal matrix elements 
vanish. 

We summarize the equations for both schemes, written in matrix form. Remember, the 
function fy defines the solution for the leading order interaction term. 

I. The flow equations by Wegner Q: 

ij = \n HA- + — iJ - — (167) 

dl [ ^ Ul ^ + dl Uij ' [ ' 

with 

Uij = exp(-ZAj) (169) 

and 

fij = u^ . (170) 

II. The similarity unitary transformation by Glazek and Wilson 0: 

f^ = 0tfbiff/la+ra *s|k, (m) 

% . = ^V(b,^k-^^] (172) 

and 
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(173) 
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Also the following transformation is used 0: 

*-^( r «t^'i«-^3 (i75) 

and 

/ij = «ij , (176) 

where for both similarity transformations 

u ij = 9(\-\A tj \) (177) 

and 

u ij + r ij = l. (178) 

Also other choices for the similarity function Uij with the step behaviour are possible ||. 

Remember, the function fy defines the solution for the leading order interaction term. The 
first order equations for H and i], written through the /-function are unique for both methods 
(I. and II.) 

<l( l ) = <U l = °)j^r^- y ( 179 ) 



1 dH^ 
" [r> Ei-Ej dl 



,W 



^\l) = —ET^-ET —^ (180) 



with the connection given in eq. (|166|) in the renormalized values, and dl — > dX implied. This 
will be exploited further for the calculations in the main text. 

B Calculation of [rj^(l), H eej ] in the ee-sector 

Here we calculate the commutator [r]^(l),H en ] in the electron-positron sector. The leading 
order generator t/ 1 ) is: 



^ (1) (0 = E / (v; m (l)e\a q + V PlP3 (l)e i r^ q ) (KM+KJt pl + d^J Pl + d.jtj 

x xt^ t i{pi,P3,-q)xs 1 S q ,- {pi - P3) , (181) 

where 

%>1P3\V ~~ _ ^P1P3 ' 9piP3 — ~~Z ' ~Ti > \*-°^) 

piP3 
A P1P3 = Pi — P3 — ijpi — Pz)~ i and the electron-photon coupling 

Hee-y = E / (#P2P4 (0 £ A'«<?' + 9 P 2P4^) £ y^-g') ( & p 4 & P2 + K^~P2 + d ~P4 b P2 + ^-P4^-p 2 ) 

xxfS]{p2,P4,-q')xs 2 V,-(P2-p 4 ) , (183) 

where 

T Pl ,p 2 , 9 = 2^ ^ a 1 - a 1 —^ 184 

q+ pi pi 
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and the tilde-fields are defined in eq. (|33|). Further we use the identities for the polarisation 
vectors and spinors 

5>M = #*. x + Xa , = S ss ,. (185) 

A 



Using the commutation relations, eq. 
M)(1\ tj i _ L I - ~* ^(^i 



), and identities eq. (|185|) we have 



fo w (Q,ff. 



ee7J 



'Vprpad- 



Ps 



P2P4 



pT-pI 



i * 



0(p1 



pi) 



PlP3 yp2P4 + 



Pl~ 



(186) 



x: ( b p3 d_ p2 d- P4 b Pl 



b + d + d b -\- b + d + d b + b + d + d b * 

U P4 U '-P1-PS U P2 ' U p 3 U '-p 1 U '-P4 U P2 ' >4 -P2 P3 U Pl, 



X (X^ r KPl 5 P3,Pl -P3)X«J (xi F KP2,P4,P2 -P4)Xs 2 ) 5p 1+ p 2 ,p3+p 4 , 

where the first two terms of the field operators contribute to the exchange channel, and the 
next two to the annihilation channel. We take into account both s- and t-channel terms to 
calculate the bound states. The : : stand for the normal ordering of the fermion operators and 
(5) is the symmetry factor. The sum over helicities S{ and the 3-dimensional integration over 
momenta pi, i — 1, ..4, according to eq. (|34"D is implied. We rewrite for both channels 



where 



fa* 



M!g\h: 



~> _,.! : K'lpi,P39-p 4 -p 2 



tot-Pa) {li 

e(~(p+-p+)) 
-frl-pt) 



t^TZw^V yl-p4,-p29 pi 



1 1-P4-P29P1,P3J 

r lpi,P39-P4-P2J 



!>■>, 



xd d Pl +p 2i p 3 -|_p 4 b p3S3 d piSi d P2S2 b PlSl 



ec-fj 



M. 



(ex) 



2ij 



M (an), is \ e(p++p+) , « 

1VI 2ij l2i|(p+ +p +) \Vpi-P29- 



e(-( P ++p+)) r 



4- " /I+ ,'J!+\" {V-P4,p 3 9p l - 



P4,P3 r l-p4,Pz9p\-p2) 

r lp\-p29-P4,P3, 



!>2 



xd d Pl + P2) p 3 +p 4 b p3S3 d piS4 dp 2S2 bp lSl 



(xfJ}(PuP3,Pi -Pa)Xn) (X? 2 rf(-p 4 , ~P2,-(Pi -Ps))x 



A' 4 , 



;i87) 



M, 



(an) 



■2ij 



(xt r !(-P4,P3, ~(P1 +P2))Xs 4 ) (%S,n(Pl> -P2,Pl +P2)Xsi 



The first term in the exchange channel with pi > p% corresponds to the light-front time 
ordering Xi < x% with the intermediate state P^ = p% + (jp\ — pz)~ + p 2 , the second term 
Pi" < P^ an d x i~ > ^3 nas the intermediate state P^ = P\ — (pi — P3)~ + pj . Both terms can 
be viewed as the retarded photon exchange. The same does hold for the annihilation channel. 

Consider only real couplings and take into account the symmetry 



V- 



P4,-P2 



-Vp 



4,P2 1 



9- 



P4,~P2 



9p 



'4,P2 • 



(189) 



Then ^3S3,p 4 s 4 | [f/ 1 ', H en ] \piSx,p 2 s 2 ) , the matrix element of the commutator between the free 
states of positronium in the exchange and annihilation channel, reads 



<[r ] W,H en ]>/5, l 



M. 



2ii („+_„+ 



P1+P2.P3+P4 



-M. 



(pT-pV 



\Vpi,p 3 9p4,P2 "T Vp4,P29pi, 



VA / 



2ii W+pT) 



(v, 



9 P 4 



Pl--p2ilP4,-P3 



' Vp4,-P39pi,~P2/ 



(190) 
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We rewrite this expression through the corresponding /-functions 



Vpi,P39p4,P2 ~T Vp4,P2 9pi, 



T lpi,-P29p4,-P3 ' r lp4,-p:i9pi, 



!>■; 



-pi 



' < ^JP1,P3\'') f f1\ | ' ( Vp4,P2\'') 



. Pl)P3 



"/p4,P2VV "T 



JP1,P3\'') 



"P4,P2 



(191) 

"Jp4,-p3W l" A J7 JP\-P2\ l ) 



^Pl,~P2 dl 



A 



P4,-P3 



c// 



with A pijP2 = Pi — P2 — {Pi ~ P2)- As we have mentioned in Appendix A this form in terms of 
the /-function is universal for all unitary transformations. We exploit further this expression 
by specifying the /-function to compare the effective interactions in different renormalization 
schemes (see Appendix C). 

We calculate the matrix elements M 2 u, eq. (188), for both channels. Here we follow the 
notations introduced in 0. 

We make use of the identities 



with s = — s and xfXs' = 8 S s'] also of 



xt^v'x* = $ ij + ise 1 - 



xf^Xs = ~V2se l s , \\>' U 



,+„■*., 



■y/2. 



se: 



(192) 



(193) 



with e* = —e s and e l s £ l s , = —S S g,. 

We use the standard light-front frame, fig. (£|), 

Pl = {XP + , XP 1 - + k ± ) , p 2 = ((1 - X)P + , (1 - X)P ± - k ± ) , 

p 3 = {x'P + , x'P 1 - + k' ± ) , p A = ((1 - x')P + , (1 - x')P x - k' ± ) , 

where P = (P + , P^) is the positronium momentum. 

Then, to calculate the matrix element M 2 u in the exchange channel, we find 
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p+ \xi i T\Pi,P^Pi -Ps)Xsi] 



A 6' 
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(X — X ) X 
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X XX 
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(195) 



and 
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im 
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T[5 S2Si +im _ {-y/2)s 2 el 2 8 S2 - Si 



where we have introduced 

T[ 

Pi 



2 (k ± -k' ± y | kj( S2 ) | k'\{s 2 ) 



Xs4 

(196) 
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(k ± -k' ± y ki( Sl ) fe'KsO 



(197) 



x — x' x 
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x' 



and 



Finaly we result 
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(198) 
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"SlSs^Szl 



'xx' x s > + S1S3 S2S4 (l-x)(l-x>) 2 S2 Jj 



Whereas in the annihilation channel we calculate 



P+ [xi.J\-p4,P3,-(Pl+P2))Xs i \ = X. 



•V, 



x' 



a + a 



a-k\ 



+ im 



^<5 aaSA + im 



~3 U S3S4 



x'(l — x' 



1 — x' x'(l — x' 

■V2)s 4 eZS Si 



-a 



S 4 "S3S4 



Xs4 

(200) 



and 



p+ [xfX(Pi,-P2,Pi +P2)Xs 1 ] 



X 



*2 



a — a 

1 — x 



im- 



x 



T % 4 5 Sl s 2 - im 



x(l — x) 



x(l — x) 
V2)si£* sl 5 slS2 , 



a 



Xs\ 



(201) 



where we have introduced 



ti = 


k'Kss) 

x' 


, k'\(s 3 ) 
l-x' 


rpi 


fci(si) 


fci(si) 



(202) 



1 — x 



X 



We finally have 



+2 itAan) 



P +A M. 



2ii 



Vs 1 s 2 Vs i s 4 -L3 -±4 ~^~ lS2 3S4 lSi xx'(l — x)(l — X' 



(203) 



-\-imv2 



0s3s 4 Sl 



Sl 



*2 



x(l — x) 



-T^-ei -S M „ A S. 



S3 



3 c Si "S3S4"S1S2 , 



x'(l — a;') 



J 4 ' fc s 4 



C Generated interaction in comparison with the result 
derived in the renormalization scheme of Glazek and 
Wilson 

We derive here the generated interaction to the order 0(e 2 ) in the renormalization scheme of 
Glazek and Wilson and compare it with the one derived in the scheme of Wegner. To this 
end we use the flow equations of type II., see Appendix A, with, for sake of simplicity, the 
/-function chosen as 

f Pm ,x = e(\-\A PtPJ \), (204) 
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where A p . p/ = Y,Pi ~HVf- This gives for the second order generated interaction in \ee >-sector 

jf y eff 

"Jpipt,X v pn 



55m 



(!) o- ee 71 



/^ / ,A<[^ ) ^r 7 ]>|ee->4 



iPf, 



dA JpiPf, 



with the solution 



Vp iPf ,x - -fpiPf.xJ d\' < [r] x , ,H e x r] >| eg> 



(205) 



(206) 



where the initial condition V PiP (X —>■ oo) = is implied. We make use of the expression for 
the matrix element of the commutator [tj, H eei ] in the exchange and annihilation channels, 
eqs. (190) and (191) in Appendix B, where now the derivative is performed with respect of the 
UV cutoff, i.e. d/dX. This results for the generated interaction in both channels 



v; 



(ex) 



-e 2 Mtl 



2ii,X + 

Pi -Pi 



+ JPiPf,^ 
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X> d f P1 ,p 3 ,X' f J \ / TOO d f P4 ,P2,X' f 

d\' JP4,P2,\ ,Uj/K . JX d\' -/P1.P3,A 



,dX' 



*Pl,P3,A 



A 
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(207) 
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/■OO ( ¥pi,—p2,X / £ 

J\ d\' JPi-Pz- 
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Pi +P2 

and with f piP2 being according to eqs. ( |176p and ( 177 ) 
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A) 



^Pl,-P2 ^P4,-P3 

This result for the generated ee-interaction was obtained by the authors of [fj. It is to be 
compared with the expression in the main text for the generated interaction, eq. (50). First, 
due to the ^-functions the interaction in eq. (208) is also free of divergencies coming from small 
energy denominators. As compared with the result of flow equations, eq. (50), each energy 
denominator in eq. (208) has its own relative weight in the scheme of similarity transformation. 
We rewrite eq. (208) as the sum of two terms, the first corresponding to the result of flow 
equations, eq. (208), and a second term representing the rest, which carries different weights 
for each energy denominator. 
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Here we have used the identity 

02 1 „ 



±slr^Mk + k 



s<* -*>(£-£ 



(210) 



where in the exchange channel 



0i 



^/pi,P3,A' 



' 3 > A f /A' 



A 0?A 



»/pi,-P2,A / f ,w 



(211) 



and 



(212) 



01 + 02 — 1 — fpipzfp4P2 • 

The term corresponding to the annihilation channel is treated along the same line. 

For completeness, we rewrite the result produced by the flow equations of Wegner, eq. (50), 
as follows 
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(213) 
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This interaction corresponds to the first term in eq. (208). In the nonrelativistic approxi- 
mation the second term with the difference of energy denominators vanishes; hence it does not 
contribute to the spectrum of positronium (see main text). In general it is an open question, 
whether the second term contributes to the physical values. This would mean the two methods 
not to be equivalent in general. 

D Fermion and photon self energy terms 

We calculate here the fermion and photon self energy terms, arising from the second order 
commutator [t/ 1 ), H eej ] . 
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I. We first derive the electron self energy terms. Making use of the expressions for the 
generator of the unitary transformation r/ 1 ) defined in eq. (^7|) and of H ee -y, eq. (|29|) , we obtain 
the following expression for the commutator in the electron self energy sector 



n\VpiP29p2Pl Vp2Pl9pip2> 



F2 — Pi 

-e(pi) 9(p l ~ p p e( P t) bp plX lx sl 



M 2lj ( Pl ,p 2 )S ij 



(214) 
(215) 



Pi ~Pi 
where 

M 2ij (pi,p2) = T t (p 1 ,p 2 ,pi -p 2 )r J (p 2 ,Pi,P2 -Pi) 

and the momentum integration over p±,p 2 is implied; 1/2 stands as the symmetry factor. The 
matrix element of the commutator between the free fermion states is 

< p X , Si | [?7 (1) , H eeJ ] \pi , Si > self energy 

r 0(t) 7) I 

= - / (Vpi P2 9p2Pi ~ Vp2pi9pi P 2) 8(pt) — + 1- M 2ii ( Pl ,p 2 ) , (216) 

J P2 Pi — P 2 



where the integration J is defined in eq. fl34|). We use the expression for the generator r\ through 
the coupling, namely 



VpiP29p2Pl Vp2Pl9piP2 
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P2 : 
P1-P2 : 

brings the integral in eq. (2TE) to the standard form of loop integration 

9(k + ) 

(Vp,p-k9p-k,p ~ V P -k,p9p,p-k) 6(p + ~ k + ) M 2ii (p,p- k) 



(218) 



(219) 



According to eq. fl5"7|), the integral J L A of the commutator \r}^~\ H en ] defines the difference 
between the energies (or energy corrections) 5p~[ x — Spi A . Making use of 



GM [VpiP2 9p2Pl Vp2Pl9p 



\P%1 



Ha -■ '■ Pi -p 2 ~ (P1-P2 

we have the following explicit expression: 

, , d 2 k x dk+ 6{k + ) , (-1) 
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(221) 



t-^p,p—k 

A 



where the solution for the ee7-coupling constant was used. Therefore the electron energy 
correction corresponding to the first diagram, cf. fig. (0), is 

**fc = e 2 /^^^V-fc + ) (222) 



2(2tt) 3 k+ "^ k 
xT l (p - k,p, —k)V(p,p - k, k) 
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x(-R) 
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where we have introduced the regulator R, defining the cutoff condition (see main text), 

i? = expj- 2 (^) 2 j (223) 

(note that A P) fc = A P]P _fc). To perform the integration over k = (k + , k x ) explicitly, choose the 
parametrization 



P + 
k 



x 



(xp + , Xp L + K L ) , 



(224) 



where p = (p + ,p ± ) is the external electron momentum. Then the terms occuring in Sp ix are 
rewritten in the form 



4^ - 4± + 2 ) k\ + 2m 2 x 2 
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(225) 



Therefore the integral for the electron energy correction corresponding to the first diagram of 
fig. (12]) takes the form 
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where 

f(x) = xm 2 - ar(l - x)p 2 . (227) 

In the last integral the principal value prescription for A as x — > was introduced (see main 
text), to regularize the IR divergencies present in the longitudinal direction. 

We thus have derived the expression for the energy correction which has been used in the 
main text. 



II. We repeat the same procedure for the photon self energy. The second order commutator 
[7/ 1 ), H eei ] gives the following expression in the photon self energy sector 
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(228) 
•TrM 2i j(px,p 2 )8 q -^ x - P2 ) , 



where M 2i j(pi,p 2 ) is defined in eq. Q215|) and the trace acts in spin space; the integration over 
the momenta q, p\ and p 2 is implied. The matrix element between the free photon states reads 
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that can be rewritten after the change of coordinates according to 
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P2 ~Pl 



-k 
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in the following way 

— / {r]k,k- q 9k-q,k ~ Vk-q,k9k,k- q ) 0(k + )d(q + - k + ) TrM 2ij {k, k - q) , (231) 

q J k 

where the symmetry 

V-Pl,-P2 = ~Vpi,P2 

9-pi,-p2 = 9pi,P2 (262) 

has been used. The integration of the commutator over I in the flow equation gives rise to 
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This means for the photon energy correction 
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where the regulator R 
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(235) 



has been introduced. Define the new set of coordinates 

(q-k)+ 



x 



r 



k = ((i- x ) q + j (l-x)q x + K ± ) 
q — k = (xq + , xq L — k l ) , 

where q = (q + , q 1 -) is the photon momentum. Then the terms present in 5q± x are 
r(k, k - q, q)T\k - q, k, -q) = {q+) 2^_ x) 2 ((to - 2 + J) ^ + = 



(236) 



A ^ = r - k--(q- k)- = -f^ + £ 



*k — q,k 



(237) 



The integral for the photon energy correction corresponding to the first diagram of fig. (]3|) takes 
the form 



q + Sq 
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with 

f(x) = m 2 - q 2 x(l - x) , (239) 

and the principal value prescription, denoted by '[ ]', introduced to regularize the IR divergen- 
cies. 

This is the form of the photon correction used in the main text. 
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Figure 1: Renormalized to the second order 0(e 2 ) light cone theory (the UV cutoff is 
A). The photon momenta are x + -ordered, from left to right. The similarity function 
/p i p / ,A = exp(— A 2 . p J\ 2 ) playes the role of regulator, where A PiPf x = T l p^ —T,pJ (the index L i L 

denotes initial and '/' final states) and A PlP2 x=Pi — p 2 — (pi — P2)~ , P~ = (p±+ rn x)/P + ■ 
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Figure 2: Electron self energy: the first diagram corresponds to the commutator term [rj^\ H eey ] 
in the electron self energy sector, next two diagrams arise from the normal ordering of instan- 
taneous interactions. 
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Figure 3: Photon self energy: the first diagram comes from the commutator [r]^, H eei \ in 
the photon self energy sector, the second one from the normal ordering of the instantaneous 
interaction. 
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Figure 4: The renormalized to the second order electron-positron interaction in the exchange 
channel; diagrams correspond to generated, instantaneous interactions and two perturbative 
photon exchanges with respect to different time ordering. 
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